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Robustness of edge states and non-Abelian excitations of topological states of matter promises 
quantum memory and quantum processing, which is naturally immune against microscopic imper¬ 
fections such as static disorder. However, topological properties will not in general protect quantum 
system from time-dependent disorder or noise. Here we take the example of a network of Kitaev 
wires with Majorana edge modes storing qubits to investigate the effects of classical noise in the 
crossover from the quasi-static to the fast fluctuation regime. We present detailed results for the Ma¬ 
jorana edge correlations, and fidelity of braiding operations for both global and local noise sources 
preserving parity symmetry, such as random chemical potentials and phase fluctuations. While in 
general noise will induce heating and dephasing, we identify examples of long-lived quantum corre¬ 
lations in presence of fast noise due to motional narrowing, where external noise drives the system 
rapidly between the topological and non-topological phases. 
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I. INTRODUCTION 

At present there is significant interest and ongoing ef¬ 
fort in realizing and detecting topological phases of quan¬ 
tum matter in the laboratory [1-9]. These efforts are 
driven by both foundational aspects of our understand¬ 
ing of quantum ordered phases in many-body systems 
beyond the Landau paradigm of local order parameters 
(see e.g. Refs. [10-12]), and in particular by the promises 
to use the intrinsic robustness of topological properties 
against imperfections in quantum information process¬ 
ing [13-17]. An example is provided by Kitaev’s quan¬ 
tum wire [18] supporting a pair of Majorana edge modes 
- Majorana fermions, which show non-Abelian exchange 
statistics under braiding [19, 20] and represent a topolog¬ 
ically protected non-local zero-energy fermion. In a wire 
network, these properties can be used to create topo¬ 
logically protected qubits and gate operations [20]. The 
quest to demonstrate Majorana fermions and their non- 
Abelian properties is presently an outstanding challenge 
in quantum physics [21-24], and is the focus of a sig¬ 
nificant effort involving systems from hybrid nano-wires 
[20, 25-37] to cold atom setups [38-46]. First evidence 
for Majorana edge modes has been reported in recent 
experiments [31-37]. 

However, while the promise of topological protection 
of quantum states from microscopic imperfections may 
hold for static disorder (see e.g. Refs. [15, 16]), re¬ 
cent theoretical studies have concluded that Majorana 
qubits and braiding can be seriously affected by coupling 
to an environment [47-51], as will be the case in any 
realistic experimental scenario. The protection of Majo¬ 
rana modes in the Kitaev wire is related to protection of 
fermion parity, and quantum correlations between Ma¬ 
jorana states will be rapidly destroyed by injection or 
removal of quasiparticles [49] . Even the coupling to a fi¬ 
nite temperature bosonic bath, which preserves particle 



Figure 1. Schematic of Majorana edge modes in noisy Kitaev 
wires and braiding operations in presence of colored Marko¬ 
vian noise, (a) A Kitaev wire in the topological phase sup¬ 
ports two Majorana edge modes 7 l /r- In realistic implemen¬ 
tations, a noisy Kitaev Hamiltonian (1) containing stochastic 
parameters X(t) = {Jj(t), Aj(t), fijft)} arises, due to a cou¬ 
pling to a classical environment (blue background). The fluc¬ 
tuations can be either global [e.g. /x(t)], or local [e.g. Hj{t)\, 
depending on the physical sources of noise in specific imple¬ 
mentations. (b) Two typical examples of colored Markovian 
noise X ft ): the colored Gaussian noise (upper panel) and the 
two-state telegraph noise (lower panel), (c) Braiding Majo¬ 
ranas (say, 71 and 72 ) on a noisy T-junction in solid-state 
setting (left panel), and in a noisy atomic wire network in the 
optical lattice setup (right panel). 


parity, is predicted to result in unavoidable losses of co¬ 
herence and errors [50]. Nevertheless, as we will show in 
this paper, it is possible to identify examples with long- 
lived quantum correlations between Majorana states in 
the presence of noise. In particular, we will be inter¬ 
ested in the effects of local and global noise representing 
a parity-preserving coupling to an environment, which we 
model as a classical stochastic process. The case of local 
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noise is representative of a two-level fluctuator in a solid 
state realization of a Kitaev wire 1 , while global fluctua¬ 
tions can result, for example, from laser light fluctuations 
in cold-atom experiments. 

Our goal is thus to study the effects of noise on Ma- 
jorana correlations and braiding operations, in a regime 
ranging from quasi-static disorder, all the way to the limit 
of fast fluctuations, i.e. where the noise correlation time 
is much shorter than the relevant system time scales. Al¬ 
though coupling to classical noise will eventually always 
lead to dephasing and heating, dephasing can be sup¬ 
pressed in the fast fluctuation limit, even when the sys¬ 
tem is driven by the noise e.g. between topological and 
non-topological phases. This effect of noise suppression 
with decreasing correlation time is familiar from atomic 
physics as motional narrowing. There, increasing the col¬ 
lision rate between atoms can result in a narrowing of the 
spectral lines [55], and in our context in an increased co¬ 
herence time of Majorana correlations. In such cases, we 
will also determine the optimal conditions for braiding 
time scales - as a trade off between the requirement of 
adiabaticity of braiding, and decoherence time scales. 

The emphasis and value of the present work is on 
exactly solvable models quantum many-body dynamics 
of Majorana correlations and braiding operations in the 
presence of colored Markovian noise sources, as exempli¬ 
fied by telegraph noise, n-state jump models, or colored 
Gaussian noise. While in a solid state context this should 
be understood as phenomenological models of noise de¬ 
scribing imperfections like local fluctuators, we note that 
noise in cold atom experiments can be engineered, as in 
recent studies of Anderson localization with (static) ran¬ 
dom optical potentials [56, 57]. Atomic realizations of 
Majorana fermions may thus serve as an ideal platform to 
study the effect of time dependent disorder in a controlled 
setting by appropriate modulation of the laser beams to 
mimic various noise sources, and can thus provide a di¬ 
rect experimental counterpart to the present theoretical 
study. 

The paper is organized as follows. In Sec. II, we briefly 
describe the model Hamiltonian of a noisy Kitaev wire. 
Then, in Sec. Ill, we develop techniques which allow non- 
perturbative solutions of the many-body quantum dy¬ 
namics for colored Markovian noise with arbitrary corre¬ 
lation time. Based on these techniques, Sec. IV presents 
a study on the Majorana edge correlations and heating 
dynamics for a global noise, which stochastically drives 
the system e.g. across the boundary between the topolog¬ 
ical and non-topological phases. The case of local noise 
is investigated in Sec. V. Thereafter, in Sec. VI, we 
study the effect of colored noise on Majorana transport, 
and discuss optimal conditions for Majorana manipula¬ 


1 Decoherence of a superconducting qubit due to coupling to a 
single two-state charge-fluctuator or to a set of them with dis¬ 
tributed switching rates and couplings to the qubit, was consid¬ 
ered, for example in Refs. [52-54]. 


tions to obtain the best fidelity at a given noise. The 
braiding dynamics on a noisy wire network is analyzed 
in Sec. VII, based on the noisy T-junction architecture 
and cold-atom setup, respectively. The paper closes with 
a summary and outlook in Sec. VIII. 

II. NOISY KITAEV WIRE 

Our goal is to study the dynamics of Majorana edge 
modes of the Kitaev wire in the presence of noise (see 
Fig. 1). The relevant Hamiltonian is 

N-l 

H[X(t)\ = ^ ^—Jj(t)aj a j+i + A j(t)ajOj+i + h.c. 

3 =1 

N 

- (!) 

3 = 1 

where aj and a j are the operator of spinless fermions 
on a finite chain of N sites. Here Jj(t) is the hopping 
amplitude on the lattice, A j(t) is the pairing parame¬ 
ter and Hj(t) is the local chemical potential. We assume 
that these parameters fluctuate in time according to a 
given noise model X(t) = {Jj(t), Aj(t), /Xj(t)}, which is 
motivated by a particular physical noise source related 
to a specific implementation. In hybrid nano-wires the 
above Hamiltonian (or its continuous version) arises from 
a combination of spin-orbit coupling of electrons in the 
presence of a magnetic field, and the coupling to an s- 
wave superconductor [25-28, 31-37]. Thus the effect of a 
two-level fluctuator, for example, can be represented by a 
fluctuating local chemical potential pbj (t) on a given lat¬ 
tice site. On the other hand, a realization of the Kitaev 
wire with cold fermionic atoms in a ID optical lattice re¬ 
sults from a laser induced coupling to a molecular Bose- 
Einstein condensate, where a molecule is dissociated into 
a pair of fermions in the wire, thus realizing the pair¬ 
ing term A j [38, 39, 51]. Frequency fluctuations of the 
laser light can be understood as fluctuations of the laser 
detuning, or as a global fluctuating chemical potential 

Before investigating the dynamics of the noisy wire 
(1), let us briefly describe as a reference the proper¬ 
ties of a noise-free Kitaev’s quantum wire [18], where 
Jj(t) = J, A j(t) = A, and /Xj(t) = fi. When |/i| < 2 J 
and A / 0, the wire is in the topological phase char¬ 
acterized by a gapped energy spectrum in the bulk and 
by a pair of robust Majorana edge modes yx, = y^ and 
7 R = 1r of the form -f L/R = £7 f L /RjCj [here we 
use the Majorana representation C 2 j-i = aj + aj and 
C 2 j = —i(cbj — aj ) of the operators] with / l/rj being ex¬ 
ponentially localized near the left (L) and right (R) edges 
with the localization length Im • The pair of Majorana 
edge modes represents a non-local fermionic zero-energy 
mode ctM = (7l+^7r)/ 2, and the long-range correlations 
—i{lLlR) between the Majorana edge states, which are 
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of interest for us, are directly related to the occupation of 
this mode, = 1 -2(0!^%). The robustness of 

the correlations is therefore related to the conservation of 
fermionic parity which distinguishes the two degenerate 
ground states with empty and occupied mode For 
|/i| > 2 J, when the wire is in the non-topological phase, 
there are no Majorana edge states and all excitations are 
gapped. 

The presence of noisy components in the parameters 
of the Hamiltonian (1) results in random “shaking” of 
the system giving rise to changes in population of the 
aM-mode (accompanied by creation of bulk excitations) 
and, therefore, to the decay of the Majorana correlations. 
Large-amplitude global noise, say, in the chemical poten¬ 
tial, can also drive the system across the quantum phase 
transition between different phases - topological and non- 
topological ones, i.e. between the cases with two Majo¬ 
rana edge modes and with none of them. Understanding 
the fate of Majorana edge modes, the correlations be¬ 
tween them, and their braiding in the noisy Kitaev wire 
(1) involves the study of many-body non-equilibrium dy¬ 
namics induced by the noise, and we describe the devel¬ 
oped techniques for treating such problems in Sec. III. 


III. QUANTUM DYNAMICS IN COLORED 
MARKOVIAN NOISE 


In dynamics of quantum systems external noise ap¬ 
pears as stochastic parameters in the Hamiltonian, as 
in Eq. (1), and the time-dependent density matrix equa¬ 
tion becomes a multiplicative stochastic differential equa¬ 
tion [58, 59]. Solving for the averaged density matrix 
in presence of noise with arbitrary correlation times, re¬ 
quires the development of non-perturbative techniques. 
For fast fluctuations, when the noise correlation time is 
much shorter than the system response time (white noise 
limit), a perturbative treatment, in form of a lowest order 
cumulant expansion, results in a master equation for the 
stochastically averaged density matrix [58, 59]. However, 
with increasing correlation time, the system response will 
be sensitive to all higher order correlation functions, and 
in the limit of infinite correlation time the static disorder 
problem is recovered (as in Anderson [60], and many- 
body localization [61]). In solving for quantum correla¬ 
tions and braiding dynamics in the noisy Kitaev wire we 
will rely on Markovian models of colored noise, which 
allow an exact solution of many-body dynamics for arbi¬ 
trary correlation time of the noise. 

At the heart of our solution of quantum dynamics for 
colored Markovian noise is the generalized master equa¬ 
tion for the marginal system density operator. In brief, 
we assume a Markov process X(t) = {A a (£)}, and we 
denote by H[X(t)\ the associated system Hamiltonian. 
Our goal is to solve the stochastic density matrix equa¬ 


tion {h = 1), dtp(t) = —i H[X(t)], p(t) , with p(t) the 


density operator, for the stochastic average (p(t)) s with 


angular brackets denoting the noise average. Defining a 
marginal density matrix 

P (X,t) = ( P mX(t)-X)) s , (2) 

we can derive the generalized master equation for the 
marginal density matrix (see Ref. [62] and App. A) 

dtp (A, t) = C (X) p (A, t)-i[H (A), p (A, t)\ , (3) 

with A now a time independent variable. Here C{ A) is 
generator of our Markovian noise model, as appears in 
the differential Chapman-Kolmogorov equation for the 
conditional density P(A, t\X',t') [58] 

d t P{ A, t\ A', t') = £(A)P(A, t\X\ t'). (4) 

Thus the operator £(A) provides a complete specification 
of our noise model, and appears as a damping operator 
in the generalized master equation (3). Solving (3) for 
{p{t)) 3 = f dXp(X,t) provides us with the desired av¬ 
erage. We refer to App. B for a discussion of various 
limits of solving the above equation: this includes the 
regime of quasi-static disorder, and the fast fluctuation 
limit (master equation limit). We emphasize, however, 
that by solving (3) we obtain solutions valid for arbitrary 
correlation time. 

The most general models for classical Markovian noises 
[58] are described by diffusion processes with continuous 
noise trajectory and by jump processes with A(£) taking 
a set of discrete values {A m } (m = 0, ...,A r ). In the 
former case, jC( A) corresponds to the Fokker-Planck op¬ 
erator, and a primary example is the colored Gaussian 
noise (Ornstein-Uhlenbeck process), see upper panel of 
Fig. 1(b). In the latter £(A) reduces to a matrix C mn 
describing the jump-rates between different A m , as in 
the case of multi-state telegraph noise, see lower panel 
of Fig. 1(b) for the two-state telegraph noise. While 
the solution for a diffusion process can in principle be 
obtained in terms of the eigenfunctions of the Fokker- 
Planck operator A) [58, 62], we pursue here a much 
more convenient strategy by discretizing it to a A r -state 
jump model with large N r and properly chosen £ mn - 
“putting the noise on a lattice”. In this case the gener¬ 
alized master equation (3) takes on the form of a set of 
coupled equations for marginal density matrices p( A m , t) 

N r 

dtp(A m , t) = E h-'mn p(X n ,t)-i[H(X m ),p(X m ,t)\, (5) 

n =0 

which - at least for low-dimensional processes A (t) - is 
not significantly of more effort to solve than the original 
non-stochastic version. 

Below we will illustrate above techniques for the 
Markovian two-state jump model [see lower panel of Fig. 
1(b)], which represents the simplest but relevant example 
capturing all essential features of noise dynamics. It also 
allows a straightforward extension to multi-state jump 
models, or to noise from several fluctuators, which for a 
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large number in the sense of the central limit theorem 
approaches colored Gaussian noise (see App. C). In ad¬ 
dition, we will specialize the above equations to the case 
of a quadratic many-body Hamiltonian, as relevant for 
the Kitaev wire. 

In two-state telegraph noise, the stochastic parameter 
X(t) switches randomly between two discrete values a 
and b with rate k. Equation (4) is now a rate equation 
for the probabilities P(a,t) and P(b,t) to be in state a 
or b respectively, 


d 

- P{a,t) - 


' —K K 

- P(a,t) ' 

— r 

- P(a,t) - 

dt 

[P(b,t)\ 


K —K . 

P(b,t)\ 

.— JU 

lP(b,t) \ 


(6) 

We have (X) a = (a + b)/2 and {X(t + r),X(t)) s = 
a 2 exp (— \t\/t c ) for the mean value and first order cor¬ 
relation function, respectively, with r c = 1/2 r the corre¬ 
lation time and the variance a 2 = (a — b) 2 / 4 (using the 
notation (X, Y) s = (XY) a - (X) a (Y) a ). 

The dynamics of a Kitaev wire driven by a single 
telegraph noise can be readily derived in the Majo- 
rana representation, where the quadratic Hamiltonian 
(1) can be recast as H[X(t)\ = (1/4) YUi Ki[X{t)]ciCi y 
with hu = —hu being an asymmetric Hamiltonian ma¬ 
trix [18]. In the Majorana basis, the key quantity cap¬ 
turing the system dynamics is the covariance matrix 
= Tr[(p(t)) s f], with f iZ = (i/2)[c*,c z ]. Thus 
from Eq. (5), we find for the marginal densities 

-^T(a, t) = — i[h(a), T(a, t)] — /vT(a, t) + ftT(6, £), 

jT(b, t) = -i[h(b), T{b, t)} - «r(&, t) + «r(a, t). (7) 

The desired stochastic average is (T(t)) s = T(a,t) + 
T(b, t). Before proceeding with the general solution, we 
find it worthwhile to briefly comment on the fast and 
slow (quasi-static disorder) limit. 

In the fast noise limit, when 1/k is the shortest time 
scale, an adiabatic elimination of the fast dynamics in 
lowest order perturbation theory reduces the above equa¬ 
tions to a master equation in Lindblad form, 

|<r> s = -i [h+, <r> j - T [h _, (r) s ]], (s) 

with h± = [h(a ) db h(b)]/2. The first term describes a 
coherent evolution with the average Hamiltonian = 
(h) s , while the second term is the noise-induced damp¬ 
ing, which scales as 1/k. Suppression of the dissipation 
by fast noises is called motional narrowing in an atomic 
physics context [55], and corresponds to a Zeno effect 
[63]. If the average Hamiltonian is in the topological 
phase and thus supports Majorana edge modes, the ini¬ 
tial Majorana correlation will show slow dephasing due to 
the Zeno effect - irrelevant if a or b per se lies in the topo¬ 
logical or non-topological regime. On the other hand, in 
the limit of slow noise, i.e. when the jump rate k 0, 
we can ignore C in Eqs. (3) and (7), C —>• 0, during 
long times t <C ft -1 , so that for such times we simply 


have to perform the (quasi-)static average of the system 
dynamics. 

Below we will solve the full marginal density equa¬ 
tions (7) for Majorana edge correlation —i{^(l1r)s = 

- Y^ii fL,ifR,i{tii)s, With f L / Rt i the Majorana wavefunc- 
tion of the initial Hamiltonian in the Majorana basis. 
Our discussion will focus on the special case of local and 
global telegraph fluctuations of the chemical potential. 
As mentioned above, such telegraph noise can represent 
a local two-level fluctuator in a solid state setup, and 
global laser frequency noise in an atomic realization of 
the Kitaev wires. We would like to stress, however, that 
even though we consider here the simplest jump-like noise 
process and only the chemical potential as a fluctuating 
parameter, our conclusions remain valid for a more gen¬ 
eral scenario with colored noise and with several fluctu¬ 
ating parameters. 

We conclude this section with the remark that the con¬ 
sidered Markovian noise models give rise to first order 
noise correlation functions, which are exponentials, or su¬ 
perposition of exponentials. We note however that non- 
Markovian noise models can often be represented as pro¬ 
jections of higher dimensional Markov processes, which 
can also be solved by our techniques. Another point to 
mention is that for non-quadratic Hamiltonians in one 
dimension, the generalized master equation (3) can be 
solved using the density matrix renormalization method 
[64, 65], and the developed techniques therefore are inter¬ 
esting in a much broader context for many-body systems 
driven by colored noise. 


IV. EDGE AND BULK DYNAMICS IN 
PRESENCE OF GLOBAL NOISE 

We begin with discussing the time evolution of the 
Majorana edge correlation — i{~iLlR)s under the global 
fluctuations in the chemical potential /i(£), when it flips 
between two values fi a and fi b (the corresponding Hamil¬ 
tonian are H a and H b , respectively) with the jump rate 
k. The statistic property of fi{t) is described by the mean 
value yu ave = (fi a + fi b )/ 2, the variance a 2 = (fi a -fi b ) 2 / 4, 
and the correlation time r c = 1/2k. We denote by 
H ave = (H[fjb(t)]) s the average Hamiltonian over the noise 
realization, which in the considered case is of the form 
of a noise-free Kitaev Hamiltonian with fi = fi ave . For 
the initial condition, we take fi( 0) = fi a with \fi a \ < 2 J 
such that Hamiltonian H a is in the topological phase, 
and assume the system is in the ground state with the 
Majorana edge-mode correlation —i^lIr) = 1- For the 
value fi b , we consider three possible scenarios (see Fig. 
2(a)): (1) \fi b \ < 2J, when the Hamiltonian H b is in 
the topological phase (fluctuations within the topological 
phase); (2) \fi b \ > 2 J but \fi ave \ < 2J, when H b is in the 
non-topological phase but the average Hamiltonian H ave 
remains in the topological phase (fluctuations between 
topological and non-topological phases but on average 
staying in the topological one); and (3) \fi b \ > 2 J and 











Figure 2. Evolution of Majorana edge correlations in presence of globally fluctuating chemical potential, (a) For a chemical 
potential /i(t) flipping between fi a and fib, there exist three jump scenarios (l)-(3) in the phase diagram corresponding to the 
value of /ib- (b)-(d) Time evolution of the Majorana edge correlation function —i{^l^r) s for fi a = 0 . 2 J and (b) fib = J (c) 
/ib — 2.1 J and (d) fib — 4 J. In each scenario, we take three typical jump rates: r = 0.1 J, r = 0.7J, and r = 10 J, and consider 
a Kitaev wire of N = 60 sites with the pairing parameter A = 0.8J. We have considered two types of initial conditions, which 
correspond to the ground state of Hamiltonian H a (solid curves) and of average Hamiltonian i7 ave in the case of fast noises 
[blue dashed curves in (b) and (c)]. In the inset of (c), the dynamics for the noise with r = 10 J (blue curve) for fib — 2.1 J is 
compared with the quench dynamics (red curve), for a quench in the chemical potential from /i = 0 . 2 J to (i = fi ave - 


|/i ave | > 2 J, when both Hamiltonians H 5 and H ave are 
non-topological (large-amplitude fluctuation when stay¬ 
ing on average in the non-topological phase). The evolu¬ 
tion of —i( 7 L 7 i?)s calculated on the basis of Eq. (7) for 
these three scenarios, is shown in Figs. 2(b) - 2(d), re¬ 
spectively, in the regimes of fast, intermediate, and slow 
jump rate r. 

We see that noise always leads to decay of Majorana 
correlations, but the decay dynamics significantly de¬ 
pends on the amplitude and rate of the noise. Non- 
surprisingly, we find the slowest decay of the Majo¬ 
rana correlations in the scenario ( 1 ) when the fluctuat¬ 
ing Hamiltonian always remains in the topological phase. 
Strikingly, the dynamics in the scenario ( 2 ) shows the 
same features, even though here we have jumps between 
topological and non-topological phases: In both scenar¬ 
ios, we observe the fastest decay in the regime with 
an intermediate jump rate (r = 0.7J, black curves in 
Figs. 2(b) and 2(c)), when r is comparable with the en¬ 
ergy gap and the band width (~ J ~ A) of the Hamilto¬ 
nian H a \ whereas, in both slow (r — 0.1 J, yellow curves) 
and the fast noise (r = 10 J, blue curves) regimes, the de¬ 
cay is much slower, and the system exhibits substantial 
Majorana correlations for much larger times (> 10J -1 ). 
In contrast to this non-monotonic dependence on the 
noise rate, the decay rate in the scenario (3) grows with 
r (see Fig. 2(d)): for both intermediate and fast fluc¬ 
tuations there are no visible Majorana correlations for 
times t > 10 J -1 , although for slow noise they survive for 
a much longer time. 

We now detail the analysis on the dynamical behavior 
of —i{~iLlR)s in scenarios (l)-(3) in the slow, fast and 
intermediate regimes of the noise, respectively. We start 
with slow fluctuations when r <C J ~ A, fi a ^)- In this 
case, the quenches between H a and H 5 - occuring at ran¬ 
dom instants - on average take place after a typical time 
/c -1 , which is much larger than all characteristic time 
scales of the system. As is explained in App. D, on the 
time scale of several inverse band-widths of the system 


(~ J -1 ) after a quench, the Majorana correlations relax 
to asymptotic values which are determined by the overlap 
of the Majorana edge-mode wave-functions for the Hamil¬ 
tonians H a and H^. Such asymptotic values then remain 
constant till the next quench occurs. For the scenarios 
(2) and (3), when is non-topological and has no edge 
modes, the overlap is zero, and already the first quench 
completely destroys the correlations. We therefore have 
—i( 7 L 7 i?) s ~ exp(— Kt) for these scenarios (yellow curves 
in Figs. 2 (c) and 2 (d)). On the other hand, for the 
scenario ( 1 ) when has the zero mode, the overlap is 
non-zero. In this case, each quench reduces the correla¬ 
tions by a factor of < 1 related to the overlap [see 
Eqs. (Dl) and (D3) in App. D], resulting in a slower 
decay of the correlation —i{^(l1r)s ~ exp[— r(1 — G^t] 
(yellow curve in Fig. 2(b)). 

In the opposite regime of fast fluctuations (r \fJLb ~ 
l^a 1 5 J rsj A), the dynamic behavior is remarkably related 
to the Zeno effect and to the quench problem. In this 
case, the evolution of —i{^l1r)s can be explained based 
on Eq. ( 8 ) for the correlation matrix, which in the con¬ 
sidered case takes the form 

^)s = -i[h+,(T) a ]-^[N,[N,(r) a ]]. (9) 

Here a = |/x& — /i a |/2 as mentioned earlier, and the ma¬ 
trices and N correspond to the Hamiltonian H ave and 
the total particle number operator in the Majorana basis. 
Following from Eq. (9), we see that: 

(i) In the limit r 00 , when the second ‘decay’-term 
can be neglected, Eq. (9) describes the dynamics of cor¬ 
relations after the quench from the initial Hamiltonian 
H a to the averaged one H ave . As a result, the asymp¬ 
totic (t -T 00 ) value of the Majorana correlation function 
-i{lLlR)s is determined again by the overlap of the wave 
functions of the Majorana edge modes (see App. D), 
but now for the Hamiltonians H a and H ave . For Hamil¬ 
tonian H ave in the non-topological phase [scenario (3)], 
there are no such modes, and after the quench —i{^(l1r)s 
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Figure 3. Heating dynamics in the bulk, (a) The absorbed energy A E{t) in unit of J as a function of time for p ave = 0 
and different ft; (b) Heating rate D s in unit of J 2 as a function of ft for different p ave . Numerical results (dotted curves) are 
compared to the predictions from Eq. (11) (solid curves). For (a)-(b), we have fixed the noise variance cr = \pb — /z a |/2 = 0.1 J, 
and have chosen A — J. 


decays to zero on the time scale of the order of the in¬ 
verse band-width of H ave . For H ave in the topological 
phase, [scenario (1) with topological H 5 and scenario (2) 
with non-topological #5], this mode exists giving rise to 
a non-zero overlap and to a finite asymptotic value of the 
correlations after the quench. 

(ii) For a large but finite ft, the second term in Eq. (9) 
adds a slow decay on top of the quench dynamics, pro¬ 
viding the asymptotic behavior of —i(7L7i?)s shown in 
Figs. 2(b) and 2(c). The short-time (t < J -1 ) behavior 
of the correlations, seen in the form of damped oscilla¬ 
tions on these figures, is sensitive to the details of the 
band structures of H a and H ave . In general, the ‘closer’ 
these Hamiltonians are, the more pronounced are the os¬ 
cillations, and the less destructive effect has the quench 
on -i(7z,7ij) s . 

(iii) Note that Eq. (9) suggests the preparation of the 
initial Majorana correlations with respect to the edge 
states of the Hamiltonian H ave , not H a , which is prac¬ 
tically also more natural for fast fluctuations. In this 
case, the first term in Eq. (9) has no destructive effects, 
and —i{^(l1r)s shows slow decay due to the second term 
(dashed blue lines in Figs. 2(b) and 2(c)). This slow 
decay of Majorana correlations in presence of fast noises 
- even for a sufficiently large fluctuation amplitude out¬ 
side the topological phase - is a direct consequence of the 
Zeno effect, which reduces the dynamics with fast fluc¬ 
tuating parameters to a weakly damped dynamics with 
the averaged Hamiltonian. 

Finally, in the intermediate regime when the fluctua¬ 
tion rate is of the order of typical energy scales in the 
system, ft ^ J ~ A, p a (p), one has optimal conditions for 
pumping excitations into the system (heating), leading 
to the fastest decay of the Majorana correlations. This 
can be seen by looking at the growth in the system en¬ 
ergy under the action of noises. To illustrate this heating 
dynamics, we consider the case when p(t) jumps between 
p a = P — 5/jl and p^ = p + Sp (such that p ave = p) with 
a rate ft, focusing on the weak noise limit Sp <C J and 
choosing A = J. We assume that initially the system is 
in the ground state of the average Hamiltonian H ave [but 


still p(t = 0) = /i a ], with energy Eq calculated with the 
corresponding initial density matrix p s ( 0), and calculate 
the system energy gain A E(t) = Tr (p s (t)H ave ) — Eq for 
different p and ft. A typical time evolution of A E(t) for 
different noise jump rates is shown in Fig. 3(a). There, 
the energy is seen to grow linearly, A E(t) = D s t, after 
a short transition time (t < J -1 ). Then after a rela¬ 
tively long time (not shown) due to the small noise am¬ 
plitude Sp < J, it saturates to the asymptotic value 
which depends on p but not on ft (infinite-temperature 
state). The heating rate D s depends on both ft and p. 
Figure 3(b) shows D s as a function of the jump rate ft 
for the values of chemical potential p = 0, 2J, and 4J, 
which correspond to the topological, critical, and non- 
topological phases, respectively. We see that D s is a non¬ 
monotonic function of ft, which is small when noise is slow 
or fast (Zeno effect), and has a pronounced maximum for 
ft « 2 — 4J. The maximum corresponds to the situation 
when the inverse noise correlation time r -1 = 2ft, which 
determines the frequency width of the noise correlation 
function, lies inside the band of bulk excitations (the ex¬ 
act position depends on the band structure). 

The linear growth of the energy for the considered 
times can be understood by calculating the energy gain 
in the second-order perturbation theory (ft = 1): 


AE(t) =J2tv\M„\ 2 / 


dr exp (i£ u r)p(r) 


, (10) 


where M v is the matrix element of the number opera¬ 
tor between the ground state and the excited state \v) 
with the energy Eq + e v . For the considered global per¬ 
turbation, | z/) simply corresponds to the state with two 
single-particle excitations with momenta k and — fc, and 

the energy e v — 2E p = 2y^(2 J cos pa + p) 2 + 4A 2 sin 2 pa 
with a being the lattice spacing. Performing the time 
derivative and assuming the time t being larger than the 
noise correlation time r c = (2ft) -1 , we obtain the follow- 
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Figure 4. The effect of local noise at different locations in 
a Kitaev wire on Majorana edge correlations, (a)-(b) Time 
evolution of —i{^l^r) s when a noise £(t) arises in the local 
chemical potential at three different sites (d — 1,3,6), with 
a jump rate (a) r — 0.7 J and (b) r — 10 J. We take £ a = 0 
and £b — 0.8J, N — 60, A — 0.4J, and p — 0.4J. For 
this parameter choice, the localization length of the Majorana 
modes is Im ~ 2.5a, with a being the lattice constant. 


ing expression for D s : 

/ oo 

dr (h(t)ii(O)) exp(2 iE p r), 

-OO 



Figure 5. (a)-(b) The effect of a large-amplitude noise £(£) 
in the local chemical potential on the dynamics of Majorana 
edge correlation —i{^l^r)s, when the local noise f(£) arises 
at site d — 3 (blue curves) and d = 6 (black curves), flipping 
between £ a = 0 and £& = 4.4J. (c)-(d) Evolution of the edge 
correlation — 2 ( 71 , 7 #} s when noise induces at random instants 
a splitting in the wire, between sites 3 and 4 (green curves) 
and between sites 6 and 7 (blue curves). In (a)-(d), we take 
two different noise jump rates, r — 0.7 J (left panels) and 
r — 10 J (right panels). Other parameters are the same as 
Fig. 4. 


where the summation is over the Brillouin zone jd E 
(—7r/a, 7r/a\. With the expression \M P \ = \u p v p \ = 
(A 2 sin 2 pa)/E 2 for the matrix element and (/i(r)/x(0)) = 
(i 2j rcr 2 exp(—2 r |r|) for the noise correlation function, we 
finally obtain (k = pa) 


D « = 4 <r 


f 


dk 


r A 2 sin 2 k 


27 r r 2 


E\ 


Eh 


( 11 ) 


for the energy-growth coefficient. The above expression 
is plotted as solid lines in Fig. 3(b) and is in very good 
agreement with numerical data. 


V. LOCAL NOISE 

We next discuss the effects of the local noise on the 
Majorana correlations — 2 ( 71 , 7 #) s , which we model by 
adding a fluctuating part to the chemical potential on 
the site d, such that Pj(t) = fi + Sj,d£{t) with the am¬ 
plitude £(£) being described by the telegraph noise: £(£) 
randomly flips between £ a and £& > 0 with the rate r , 
and £(0) = £ a . In this case, an exponential localization of 
the Majorana modes near the edges with the localization 
length Im leads to a very strong dependence of the effects 
of the noise on d. This is because the decay of the Ma¬ 
jorana correlations is caused by noise-induced changes in 


the population of the associated non-local fermionic zero- 
energy mode, with the corresponding matrix element be¬ 
ing proportional to the value of the edge-mode wave func¬ 
tion on the noisy site d. The results of our calculations 
for intermediate and fast noises, Figs. 4(a) and 4 (b), re¬ 
spectively, clearly show this dependence: the decay is the 
fastest for d < Im, for d ~ l M it is already substantially 
less, and is exponentially small for d > Im- 

Similar to the case of a global noise, the reduction of 
the decay in the fast-noise regime (r J, A, p) is due 
to the Zeno effect, as follows from Eq. (8) which now 
takes the form 

± { T) s = -i[h + ,(T) s ]-^[n d ,[n d ,(T)s}], 

where a = |£& — £ a |/2, the matrix rid corresponds to the 
local density operator a\ad at site d in the Majorana ba¬ 
sis, and the second term adds a slow decay on top of the 
result of the quench described by the first term. The 
average Hamiltonian in this case contains the static im¬ 
purity potential Vd = (£ a + £b)/2 on the site d, which 
results in just modification of the Majorana edge modes, 
and hence in a finite asymptotic value of —i{^(l1r)s after 
the quench. The small (~ ft -1 ) decay rate is extra re¬ 
duced, as compared to the global noise, for d > Im due to 
the smallness of the edge-mode wave functions on the site 
d. For d > Im, this gives an exponentially small decay 





























rate such that the Major ana correlations are practically 
immune to the noise (Fig. 4). 

If we increase the amplitude of the noise to larger val¬ 
ues, see Figs. 5(a) and 5(b) for £5 = 4.4J, the oscillating 
behavior for fast noise become much more pronounced. 
In this case, the strong static impurity potential in the 
average Hamiltonian splits the wire into two pieces which 
are weakly coupled through the impurity site (Josephson 
junction). If the coupling would be zero, one would have 
an extra pair of Majorana modes at the edges adjacent to 
the impurity, and the corresponding fermionic zero-mode. 
For a small but finite coupling, the energy of this mode 
is finite, and the oscillations seen in this case correspond 
to the energy of this mode. Similar oscillatory behavior 
is observed when the fast noise randomly splits the wire 
into two parts, for example, when the local hopping Jd 
and the pairing A^ amplitudes (between site d and d+ 1) 
jump simultaneously between the finite values and zero, 
as shown in Figs. 5(c) and 5(d). In this case of fast noise, 
the average Hamiltonian has a “weak link”, and the os¬ 
cillation frequency seen in Fig. 5(d) corresponds to the 
energy of the fermionic mode localized at this link. Note 
that the amplitude of the oscillations is related to the 
overlap between the Majorana edge mode and the wave 
function of the low-energy fermionic mode localized on 
the “defected” site or link, and rapidly decreases with 
increasing distance between the modes. 

Note that for the ideal Kitaev chain (J = A and 
fi m 0), the two Majorana modes jl = c\ (7# = 
C 2 n) locate on the leftmost (rightmost) sites, such that 
~^{1l1r)s = — (r s )i,2iv- In this case, the dynamics of 
the Majorana correlations is completely uncoupled from 
that of the bulk - the Majorana correlations in this ideal 
case are absolutely insensitive to what happens in the 
bulk. 


VI. COMPETITION BETWEEN NOISE AND 
ADIABATICITY IN MAJORANA TRANSPORT 

Equipped with above understanding of the non¬ 
equilibrium dynamics of a noisy Kitaev wire, we now dis¬ 
cuss the effect of a local noise on the Majorana edge cor¬ 
relations during the adiabatic transport [20, 29, 30] - an 
essential building block for the braiding operations. As 
we will see, in the presence of a noise, the adiabaticity of 
the transport - required for preserving the information 
encoded in the Majorana correlations - confronts with 
the finite life-time of the correlations, and the competi¬ 
tion of these two factors establishes an optimal operation 
time. 

Following Ref. [20], we will move the left Majorana 
edge mode by “pushing” it to the right via adiabatically 
switching on local potentials on the corresponding sites. 
For example (see Fig. 6(a)), the move of 7 l from site 1 
to site 2 can be achieved by applying the local potential 
\{t)V at site 1 [an extra term \{t)V a\a\ in the Hamil¬ 
tonian], where V 2 J and \{t) increases monotonically 



Jt JTf 


Figure 6 . Fate of a Majorana moving through a local noise, 
(a) A move of the Majorana edge mode 7 l from site 1 to site 2 
by adiabatically tuning the local potential V in a time interval 
Tf. (b)-(e) Time evolution of the correlation between the edge 
Majorana operators — i(c 2 j-iC 2 N) s for j = 1,2, 3,4, when 7 l 
is adiabatically transported from site 1 to site 4 through a 
local noise £(£) at site d — 2 with (b) ft = 0 (c) ft = 0.5 J 
and (d) ft = 10J. (e) The time evolution of the correlations 
between the Majorana edge modes —i{^L{t)^(R)s for different 
ft. ( JTf = 25 for (b)-(e)). (f) The remaining Majorana corre¬ 
lations — 2 ( 7 ^ 7 ^) f as a function of Tf for different ft, after 7 l 
is moved to site 4 (in total times of 3 Tf). For (b)-(f), we have 
chosen £ a = 0, £& = 0.8J, N — 40, A = 0.8J, and r — 0.2J. 


from A (t < 0) = 0 to A (t > Tf) = 1 during the time in¬ 
terval [0, Tf] with Tf being much larger than the inverse 
energy gap, Tf J -1 , A -1 , /x _1 . (In our calculations 
we use A (t) = sin 2 [(7r/2Xy)t].) Further moves can be 
achieved by applying the same protocol successively to 

sites 2, 3,-In Fig. 6(b) we show the evolution of the 

correlations —i(c 2 j-iC 2 n) with j = 1,2, 3,4, during the 
adiabatic move of 7 ^ from site 1 to site 4 for Tf = 25 J -1 
(in total time 3 Tf) and in the absence of the noise. (The 
correlation between the actual edge modes 

remains unchanged from its initial value 1, see also the 
curve with ft = 0 in Fig. 6(e)). 

The same correlations in the presence of the local noise 
Vd[£(t)] = £(t)rid at site d [here £( t ) flips again between 
£ a and £5 at a rate ft] are shown in Figs. 6(c) (ft = 0.5 J) 
and 6(d) (ft = 10J) for d = 2. The corresponding be¬ 
havior of the correlation —i{ r )L{t)lR)s between the edge 
modes [7 l(T) is moving and 7 r is fixed] is presented in 
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Figure 7. Majorana braiding on a noisy T-junction, where a local noise £( 2 ) occurs in the common joint connecting three wire 
segments. Upper panel: schematics of an exchange of two Majoranas 71 and 72 following Ref. [20]: (i) Initially the horizontal 
wire is topological supporting two Majorana edge modes 71 and 72 at the ends, while the vertical wire is non-topological. (ii) 

71 is moved through the junction to the bottom of the vertical wire, while 72 is fixed at the right end of the horizontal wire, 
(iii) 72 is moved all the way to the left end of the horizontal wire, while 71 is fixed, (iv) 71 is moved upward through the 
junction and then rightward, (v) At the end of the exchange, the wire returns to its original configuration, with 71 —72 and 

72 —> — 71 . Lower panel: Evolution of the Majorana edge correlations — 2(71 ( 2)72 (2)) s in each step of the exchange processes, 
when a local noise £(2) in the junction jumps between 0 and 0.7J at a jump rate k = 0.6J and k — 10 J, respectively. For other 
parameters, the pairing parameter of the horizontal and vertical wires are chosen as = iA y — 0.7J. We take p — 0.1J in 
the topological segment and p — —4 J in the nontopological segment of the wires, and JTf — 18. 


Fig. 6 (e). We clearly see deviations from the noise-free 
case, which are significant for k = 0.5J and small for 
k = 10J (Zeno effect), but these deviations take place 
mostly when the Majorana mode moves to the noisy site 
(for 2 from 0 till Tf in the considered example). After 
that, the Majorana correlations do not exhibit any vis¬ 
ible decay and repeat the pattern of the noise-free case 
but with the reduced amplitudes. This behavior follows 
from the localized character of wave functions of the Ma¬ 
jorana edge states: the correlations between them are 
influenced by the local noise only when the moving Ma¬ 
jorana mode and the noisy site are within the localization 
length Im- (Note that the extend of the edge-mode wave 
function in the ‘non-topological’ part of the wire - the 
sites with non-zero local potential V is also non-zero 
but very small for V 2 J. As a result, a noisy site in 
this part of the wire has no effect of the correlations.) 

Above results also imply that, in order to minimize 
the destructive effect of the noise on Majorana edge cor¬ 
relations, the move through the noisy site has to be per¬ 
formed with the fastest speed - the requirement which is 
opposite to the adiabaticity condition for the transport. 
As a result, there exists an optimum speed of transport 
(optimum Tf) for each k>. This is illustrated in Fig. 6 (f) 
which shows the remaining correlations (after the total 
move) as a function of Tf for different k: the decrease 
in the correlations for small Tf is due to non-adiabatic 
effects, while for large Tf it is due to accumulated ac¬ 
tion of the noise. The proper choice of Tf can substan¬ 
tially reduces the loss of correlations, especially in the 
intermediate-noise regime. 


Similar consideration is also applicable to the global 
noise. However, in this case the destructive effect of 
the noise is independent of the position of the Majorana 
modes, so that an entire time of the operation should be 
within the life-time of correlations, see Fig. 2. 


VII. FIDELITY OF BRAIDING IN A NOISY 
WIRE NETWORK 

Finally, we study the effects of the noise on the Majo¬ 
rana braiding (exchange) - the operation which for the 
two modes 71 and 72 corresponds (up to a phase) to the 
unitary operator U = exp [—(^ 4 ) 7172 ] and results in the 
transformation 71 -A 72 , 72 —>• — 71 , showing non-Abelian 
character of Majorana fermions [19, 20]. We consider two 
proposed braiding scenarios: (i) in the T-junction for the 
solid-state heterostructures, see Ref. [20], and (ii) in the 
wire-networks for cold-atom systems, see Refs. [44, 45]. 

We first consider the T-junction with two Majorana 
edge modes 71 = 7l and 72 = 7 r which we braid by 
moving them, see Fig. 7, in accordance with the pro¬ 
tocol from Ref. [20]. We choose initially — 2 ( 7172 ) = 1 
and follow the evolution of this correlation during the 
protocol. Without noise, it remains unchanged during 
the entire braiding, provided we move the modes, say 71 , 
adiabatically (we choose Tf = 18J -1 ). In the presence of 
local noise, from the previous results we expect the de¬ 
crease of — 2 (71 ( 2)72 ( 2 )) s each time when the Majorana 
mode passes the noise site. This is demonstrated in Fig. 
7 for the case when the noise-source is located at the 
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Figure 8 . Braiding dynamics of Majoranas in a noisy atomic 
wire network, (a) Cold-atom implementation of braiding 
based on Refs. [44, 45], where the Majorana modes 7 ^ and 
on the left ends of the two wires are exchanged, (b)-(e) 
Time evolution of the Majorana correlation functions for (b) 
noise-free case (c) a local noise at site S 3 of the upper wire 
with k — 0.5 J, and (d)-(e) a local noise at site s\ of the 
upper wire with (d) k = 0.5J and (e) k = 10J. For (c)- 
(e), we choose £ a = 0 and £& = 0.6J. For other parameters, 
A (C = A (0 = 0 . 7 J and /x (n) = // z) = 0 , A (n) = N (z) = 32. 
Following Refs. [44, 45], the braiding is realized in four 
steps, where we choose the operation time tf for each step 
as Jtf = 25 for (b)-(e)). 


common point of the three legs forming the T-junction: 
here the Majorana modes have to cross the noisy site 
three times, and each crossing results in the decrease of 
For the noisy site located in one of the 
legs, only two crossings occur with the two corresponding 
drops in —i( 71 ^) 72 ( 0 )s? resulting in a higher fidelity of 
the braiding operation. 

For the braiding in an atomic wire network, we con¬ 
sider two wires (see Fig. 8 (a)): the upper one (u) 
and the lower one (1), each having a pair of Majorana 
modes ( 7 ^, 7 ^) and ( 7 ^, 7 ^). The braiding protocol 
from Refs. [44, 45] involves operations only on one side 
(say, left) of the network, and the mode to be braided, 

71 = 7 ® and 72 = 7 ^\ are also located on the same 
side. As a result, the protocol will be only sensitive to 
noise located close to the left side of the network. Fig. 
8 (b) shows the evolution of correlations between various 
Majorana modes during braiding in the absence of the 
noise. The evolution of the same correlations with the 
telegraphic-noise source with k = 0.5 J on the third site 
(S 3 ) and on the first site (s\) of the upper wire are pre¬ 
sented in Figs. 8 (c) and 8 (d), respectively. They clearly 


show the above mentioned feature of the protocol. No¬ 
tably, for a fast noise, even when the noise source is lo¬ 
cated on the first site s \, one has much less noise-induced 
decoherence whence higher fidelity, see Fig. 8(e). 


VIII. CONCLUSIONS AND OUTLOOK 


To summarize, we have studied the decoherence of Ma¬ 
jorana edge correlations and braiding dynamics in colored 
Markovian noises preserving parity symmetry. Our anal¬ 
ysis relies on a technique for solving quantum many-body 
dynamics when the system parameters undergo local or 
global fluctuations modelled by classical stochastic pro¬ 
cesses with arbitrary correlation time. Our studies on 
noisy Kitaev wires show that, while the noise always gives 
rise to the decay of the correlations between Majorana 
edge-states, there are several parameter regimes where 
the life-time of the correlations remains sufficient for 
quantum manipulations with Majorana fermions, even 
without error-corrections. This includes the cases of slow 
global noise and generic local noise in the bulk, and in 
particular, the case of fast noise where decoherence can 
be suppressed due to motional narrowing, also known as 
the Zeno effect. These results further allow us to opti¬ 
mize the manipulation protocols of Majoranas in both 
the solid-state and cold-atom settings. Our presentation 
is for two-level telegraph noises in chemical potentials, 
but the essential features of noise dynamics are also seen 
in the colored Gaussian noises (the lattice model), and 
in other types of noises, e.g. the phase fluctuation in the 
pairing parameters. 

The present study and the development of techniques 
to treat the effect of noise from static disorder to the 
rapid fluctuation limit, should also be seen in the broader 
context of dynamics of correlations in an interacting 
many-body quantum system in the presence of random 
fluctuations. The effect of random fluctuations, either 
as spatial disorder or temporal noise, on the proper¬ 
ties of quantum many-body systems is a long-standing 
and important problem. Spatial static disorder under¬ 
lies phenomena like Anderson localization, and quantum 
many-body localization-delocalization transition in the 
presence of interaction, with typically short-range cor¬ 
relations in the localized phase. The time-dependent 
random fluctuations introduce temporal decoherence and 
possible heating, resulting in finite temporal correlations. 
Whether the combination of random fluctuations and 
interparticle interactions could lead to some interesting 
long-range dynamics in the space-time domain, is an open 
and intriguing question. The considered model, being 
formally quadratic, implicitly contains the effects of inter¬ 
particle interactions in the form of a paring term which is 
responsible for the existence of the non-trivial topological 
phase with non-Abelian Majorana states. Due to topo¬ 
logically protection, these states and the long-range cor¬ 
relations between them survive the static disorder, and, 
therefore, the considered model provides a simple and 
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tractable example from a very special class of topological 
system, both interacting and non-interacting, with corre¬ 
lations robust against static disorder. The results of our 
paper provides therefore a possible scenario for behavior 
of such systems in the presence of a temporal noise. 
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trix p(Y, £) in the main text. Denoting A[X(t)\p(t) = 


—i 


H[X(t)\,p(t ) , 


for each noise realization we can 


solve the multiplicative stochastic equation p(t) = 
A[X(t)\p(t) for the density matrix p(t) given an ini¬ 
tial one p(to) at to with X(to) = Xq. The formal 
solution can be written as p(t) = Pn(t), with 

Pn(t) = / t * dt n ... / t * 2 dti n" =1 A[X(ti)]p(t 0 ) for t > t n > 
... > to and po(t) = p(to). Thus from Eq. (2) we have 
p(X,t) = En=OPn(X,t) With p n (X,t) = (p n (t)5(X — 
X(t))) s , where the stochastic average can be straight¬ 
forwardly performed using joint probability densities 
P(Y, t; Y n , t n ;...; Yo, to). From the defining property of 
a Markov process, that is the factorization property of 
the conditional probability densities [58] , we can write 


Appendix A: Derivation of the generalized master 
equation 

Following Ref. [62], here we derive the general¬ 
ized master equation (3) for the marginal density ma- 


pn(X,t) = f dt n .. f dti f dX n .. [dX 0 A(X n )...A(X 1 )p(t 0 )P(X,t\X n ,t n )P(X n ,t n ; ..;X 0 ,t 0 ). 
Jto Jto J J 


(Al) 


Here P(Y, t|Y n ,t n ) is the conditional probability for 
finding X(t) = X given X(t n ) = X n at the earlier 
time t n < t. The significance of Eq. (Al) is that 
the time dependence in the stochastic parameter X(t) 
is now transferred into P(Y, £| Y n , £ n ). Then by using 
the Chapman-Kolmogorov equation (4) for the evolution 
of P(Y,*|Y',*') (t' < t), together with P(X,t\X',t) = 
5{X — Y'), we obtain 

Pn(Y, t) = A(Y)p n _ 1 (Y, t) + £(Y)p n (Y, t), (A2) 
and p(Y,t) = J2n=o Pn{X,t) readily gives Eq. (3). 

Appendix B: Fast fluctuation and quasi-static limits 

Below we solve the generalized master equation (3) 
for the average density matrix p s (t) = (p(t)) s in two 
limiting cases of a stationary colored Markovian noise: 
(1) the fast fluctuation limit and (2) the quasi-static 
limit. (As above we will use the notation A(Y)p(Y, t) = 
-i[H(X),p(X,t)].) 

Fast fluctuation limit- In this case, a master equation 
for p 8 (t) can be derived using the eigenfunction expan¬ 
sion method [58, 62]. Denoting the left (right) eigenfunc¬ 
tions of the noise operator C as P\(X) [Qa(Y)], with 
CPx(X) = -XP x (X) [£t Qx(X) = -XQx(X)], we ex¬ 


pand 

p{X,t) = Y J Px{X)C x {t). (Bl) 

A 

Note that Po(X) = P s (X) represents a stationary dis¬ 
tribution with jCP s ( X) = 0, and Qo(X) = 1. De¬ 
noting (Qaj AP\r) = f dXQ x {X)A{X)P x ,{X) and us¬ 
ing the orthogonal condition (Q\,P\r) = S X) x', the ex¬ 
pansion coefficient C\(t) in Eq. (Bl) is derived as C\ = 
J dXQx{X)p(X,t). Importantly, we identify 

Ps(t) = C 0 (t) = J dXp(X,t), (B2) 

which is just the desired average density matrix. 

We thus want to derive p s = Co(t). Substituting Eq. 
(Bl) into the generalized master equation (3), we find 

Ps = A a ve p s + ^^(Qo, APx)Cx, (B3) 

A/0 

C A = -A Cx + ( Qx,AP 0 ) Ps + (Qx,APy)Cy.(B4) 

A'#0 

with A ave = J dXA(X)P s (X). For fast fluctuations 
when the damping rate ^ A is large, we can eliminate 
the fast dynamics of C x (t) [see Eq. (B4)] on a time scale 
t A -1 using the technique of adiabatic eliminations 
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[58]. In doing so, Equation (B3) becomes (in linear order 
of 1 /A) 

p s (t) = (A ave + V) p s (t) . (B5) 

Here the operator V is defined by Vp s = 
fo dT(A(r), A(0)) a p s , with (A(r), A(0)) s = 

Sa/oLT dXA(X)P\(X)] 2 e~ Xr the stationary vari¬ 
ance [58] and t A -1 . Equation (B5) is the familiar 
master equation: the first term corresponds to a coherent 
evolution with the average Hamiltonian H ave ; the second 
term describes a damping dynamics with the operator V 
determined only by the second order noise correlations. 
As an illustration, consider H[X(t)\ = H ave + X(t)H\ 
[thus A[X(t)\ = A ave + X(t)Ai]. For the exam¬ 

ple of noises X(t) with exponential correlations, 
(X) s = 0 and (X(t),X(0)) s == cr 2 exp(—|r|/r c ), we have 
V = g 2 t c A\A\ for t r c when Eq. (B5) becomes 

ps = —i[H ave , p s ] - a 2 r c [Hi, [Hi,p s ]]. 

Similar equation arises in the main text in the concrete 
example of a fast two-state telegraph noise [see Eq. ( 8 )]. 

Quasi-static limit- We now turn to the quasi-static 
case, when the relevant times (say, time for experimental 
observation) are much shorter than the noise correlation 
time, t <C r c . For this time the noise distribution is 
effectively frozen to the initial one, and hence we ignore 
C in Eq. (3) when it reduces to 


The stochastic property of the above noise X(t) is 
described by a probability density P(X m ,£) for finding 
X(t) = X m at time £, i.e. 

(C2) 

Here P(a'/b', t ) is the probability for a telegraph being in 
state a'(b') at time t with its evolution given by Eq. (7). 
Hence we obtain 

N r 

dtP(X m ,t) = E d^mn P(X n ,t), (C3) 

n =0 

with — tvN r , £ m?m+ i — ft ( 171 T 1), Em,m—1 — 

ft (N r — m + 1), and « 0 for n ^ m, m d= 1. 

For large N r , the distribution (C2) approaches a Gaus¬ 
sian distribution as ensured by the central limit theorem 
[58], and Eq. (C3) represents the Fokker-Planck equa¬ 
tion for an Ornstein-Uhlenbeck process [58]. To see this, 
we note that X m+ i — X m = b' — a' ~ 1/y r X r (with 
fixed a of the noise), so that when N r oo we can re¬ 
place X m with the continuous variable X and expand 
P{Xm± l) in terms of AX = b' — a' as P(X ± AX) « 

P(X, t) ± (6' - a')d x P(X, t) + { -^p^ d 2 x P(X, t). In view 
of m = ( N r b f — X)/(b r — a') from Eq. (Cl), we obtain 
from Eq. (C3) that 


p(X, t) = -i[H(X),p(X, t)} = A(X)p(X, t). (B 6 ) 

Given an initial p(X, 0), solution of Eq. (B 6 ) gives the 
average density matrix for times t < r c as p s (t) = 
f dX exp[A(X)t]p(X, 0) (valid to lowest order of r c _1 ). 


Appendix C: A lattice model for colored Gaussian 
noise 

Here we present a lattice model for colored Gaussian 
noise X(t) (Ornstein-Uhlenbeck process [58]), character¬ 
ized by a mean value (X(t)) s and a variance (X(t + 
r),X(t)) s = a 2 exp(—|r|/r c ). The basic idea is to form 
a multistate noise with N r independent two-state tele¬ 
graph noises: X(t) = X r (t). Each telegraph X r (t) 

flips between a' and b' at a rate ft, with (X r ) s = (a'+b')/ 2 
and (X r (t + r), X r > (t)) s = 5 r ycr' 2 exp(—^r|/r c ). Here 
a' 2 = ( a' — b') 2 /4 and r c = 1 /2ft as in the main text. 
Thus by construction we have (X(t)) 3 = X r (X r ) s , and 
a 2 = N r a' 2 . The noise X(t) can be viewed as result¬ 
ing from many independent two-level fluctuators with the 
same jump rate, so that the instantaneous value of X(t) 
switches randomly among N r + 1 discrete values {X m } 
(m = 0,.., N r ) with 

X m = mo! + (N r - m)b'. (Cl) 

We remark that, for given (X(t)) s and a of the noise 
X(£), the values a' and b' of each telegraph are deter¬ 
mined from scaling relations: a! = (X(t)) s /N r — a/^/Nr 
and b' = (X(t)) s /N r + a/^/W r . 


dtP(X,t) = 


2 «£c (x - + \ D & 


P(X,t). 


This is the Fokker-Planck equation governing the 
Ornstein-Uhlenbeck process [58], where we identify a 
drift velocity 2ft, and a diffusion constant D = 4ftcr 2 . 

Thus the quantum dynamics of a system in (dis¬ 
cretized) colored Gaussian noise is governed by the gen¬ 
eralized master equation of a form (5) with X m and jC mn 
given by Eqs. (Cl) and (C3), respectively - for large 
X r , it conveniently approximates Eq. (3) for a colored 
Gaussian noise when C(X) is the Fokker-Planck operator 
corresponding to an Ornstein-Uhlenbeck process. 


Appendix D: Asymptotic long-time Majorana edge 
correlations in a quenched Kitaev Chain 


Here we derive the asymptotic Majorana edge-mode 
correlation at large times (t Too) after a global quench 
in the chemical potential of a Kitaev Hamiltonian, with 
the quench from po to P/- Specifically, we assume the 
system is initially in the ground state |0) of Hamiltonian 
H(p o) with |po| < 2 J in the topological phase, support¬ 
ing two Majorana edge modes 7 l/r. Suppose at time 
t = 0, the Hamiltonian is globally quenched from H(p 0 ) 
to H(pbf) with \pbf\ < 2 J. Our goal is to derive the 
asymptotic (t —» + 00 ) Majorana edge correlation 


G oo = lim ( 0 

t —^--(-oo 


[~i lLlR \ e -«*(w)t| 0 y (Dl) 
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Figure 9. Long-time asymptotic Majorana edge-mode cor¬ 
relation Goo (Dl) as a function of the post-quench chemical 
potential (i /, when a Kitaev Hamiltonian is globally quenched 
from H(tio) to H(/if), with |/io|, \(if\ < 2 J. Blue line - nu¬ 
merics, red line - Eq. (D3). We take (i o = 0.3J, N = 134, 
A = 0.72 J. 


In the Majorana basis, Eq. (Dl) can be written as 
Coo = Eijli with = lim t ^ +00 -i(cm), 


and f^j R l describes left (right) Majorana modes of the 
initial Hamiltonian. As an illustration, we calculate G^ 
for i = 2ji — 1 and l = 2j 2 with j± ~ 1 (near the left edge) 
and j 2 ^ N (near the right edge) such that c 2 j 1 -i = aj ± + 
and c 2j2 = —i(a j2 — atj. Using Bogoliubov trans¬ 
formation dj = we diagonalize the 

post-quench Hamiltonian H(fif) as H(fif) = ^2 
(up to unimportant constant) with being the quasi¬ 
particle operators for the Hamiltonian H(/if). With 
this we obtain e lH ^^ t aje~ %H ^ tI, ^ t = + 

, and similarly for After 

substituting this into the expression for G ^, and neglect¬ 
ing oscillating terms (the ones containing /^/^, 
and with (i ^ u) which are averaged to zero after 

times larger than inverse band-width, we obtain 


^2.71-1,2,72 ~ y! + -®jij 2 /x(0|^/x^|0)) 5 


with = (ii* iAt + i;* lM )^ 2 M + (^y/x + and 

Bjijzn = (^ji/x + v jin) v j 2 n (Vji/j' + u ni±) v hv-' To ca ^“ 
culate the correlations in the above expression, we use 

the second Bogoliubov transformation aj = ^ a u^°J a a + 

u j^* a t which diagonalizes the initial Hamiltonian as 

H(/i o) = Xlcr e cr a )j a (j (again up to unimportant con¬ 
stant) in terms of quasiparticles operators a a (<aj.), such 
that <a a |0) = 0. Then, by virtue of the relation = 


+ D* a al, with C„ a = + v^vfj 

and = Ej + Vj^ufJ, we arrive at 


/OOO 

1,2.72 


{Dual I^WI ) • 

H a <j 

(D2) 


Equation (D2) contains contributions from both the 
edge (/i = M ) and bulk modes (fi ^ M ) of H(/if) 
(|/ij | < 2J). Due to the gapped energy spectrum, the 
bulk contribution (modes with /i ^ M ) to the corre¬ 
lations between the edges are exponentially suppressed, 
and therefore can be ignored in Eq. (D2) in the thermo¬ 
dynamic limit. On the other hand, for the edge con¬ 
tribution (/i = M), we use Majorana wave functions 
9L/R,j °f Hamiltonian H(fif) in the fermionic represen¬ 
tation: gLj = UjM + VjM and qrj = UjM — VjM- 
Keeping in mind the localization character of the Ma- 
joranas at edges, we have Aj 1 j 2 M ~ gLj 1 9Rj 2 and 
1 j 2 M ~ —9L,j 1 gR,j 2 - Equation (D2) is then simpli- 

fied as Gf h _ 12j2 « gL,j 1 9R,j 2 H (r (|CW| 2 ~ l-CWl 2 )- 
This expression further involves contributions from the 
edge a = M and bulk modes a / M of initial Hamil¬ 
tonian H(/jLq). For same reasons discussed earlier, we 
ignore the exponential small bulk contributions (a ^ 
M). For the rest edge contribution (cr = M), we write 

u fl “ (1/2)(5l3 +9rj) and = (1/2 - g ( °\) 
with g^/R j the initial Majorana wave functions of Hamil¬ 
tonian H((i 0 ) in the fermionic basis. Thus by us- 
ing C M m = \ (9lj9lj + 9rj9rj) and D M m = 

5 Ej (9L,j9 ( Lj - 9rj9rj), we find 


^2/1-1,2j 2 


9L,j 1 9R,j 2 


' 3=1 


][ 


N 

E ( 0 ) 

9RJ29rj 2 

3 = 1 


Calculations of other components of G^ are similar, 
and we finally get 


Gqq — 


= 


3=1 


N 

3=1 


1 2 


(D3) 


We see that, after the quench, the Majorana edge corre¬ 
lation in the long time approaches an asymptotic value, 
which is determined only by the overlap between the 
wavefunctions of the Majorana edge modes of H(/jLq) and 
(i.e. g^/Rj and gL/R,j )• Figure 9 shows numerical 
results for Goo (blue curve) as a function of /if , which are 
compared to predictions from Eq. (D3) (red curve). A 
good agreement is clearly found, with deviations only ap¬ 
pear near the critical point when the energy gap becomes 
small. We thus conclude that Majorana edge correlations 
relax to a finite value after a quench within the topolog¬ 
ical phase, which decreases with /if in the post-quench 
Hamiltonian and eventually varnishes for (if =2 J. 
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